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California Institute of Technology 
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{0_>ATP, vi}, 
{ADP^0, v2}, 

{PFK + ADP 2 *± PFKA, k3, k*3}} 
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• Phosphatases act on kinases in scaffold XFirst and second phosphorylation rates 
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gms = interpret [[ 

{0mC, Yi, kd}, 

[0^0, nill[vmax Yd, khalf kd]}, 
(Comp[H] h> H, hill[Ymax -» Yi, khalf ->• Ki]}, 

hill[Ymax -» Y2, knalf-*K2]}, 
{Comp[X]h+X, hill[Ymax YH3 *M[t], khalf -»K3]} 
{Xh0, liill[Ymax -» Y4, knalf-*K4]} 
}] /. {Y1-* Yni*C[t] /(Kc + C[t])>; 
Print [gms]; 

run[gis, timeSpan-* 500, MaxStep3 5000, 
initialConditions {C[0] == .1, X[0] == . 1}, 
plotYarianles -> All, 
rules -* { 

Ki-» 0.001, K2-> 0.001, K3-»100, K4 -» 100 , 

▼i-» 0.025, lc-» 0.3, kd-» 0.01, YHi -»• 3, 

Y2-* 1.5, YH3-»i, Y4-+ 0.5, Yd -> 0.25, Kd 0.02 

checkConstants False 

] 
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FIG. 10 
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Receive initial condition values, process parameters, and a user 
representation of the biological network, wherein the user representation 
is input using one or more of a series of biological network canonical 
input forms, wherein each canonical input form is based on a type of 
biological process in the biological network. 



Generate a series of mathematical equations in an 
equation output canonical form based on the input 
representation of the biological network and the 
process parameters. 
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Numerically solve the series of mathematical equations using the initial 
condition values and the process parameters, to generate a value or a table 
of values as a function of time for one or more output functions of the 
biological network, thereby simulating the biological network. 
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